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ABSTRACT 

We investigate the evolution of the surface properties of models for rotating massive stars, i.e., 
their luminosities, effective temperatures, surface rotational velocities, and surface abundances 
of all isotopes, from the zero age main sequence to the supernova stage. Our results are based on 
the grid of stellar models by Heger et al. (2000), which covers solar metallicity stars in the initial 
mass range 8 — 25 M©. Results are parameterized by initial mass, initial rotational velocity and 
major uncertainties in the treatment of the rotational mixing inside massive stars. 

Rotationally induced mixing processes widen the main sequence and increase the core hydro- 
gen burning lifetime, similar to the effects of convective overshooting. It can also significantly 
increase the luminosity during and after core hydrogen burning, and strongly affects the evo- 
lution of the effective temperature. Our models predict surface rotational velocities for various 
evolutionary stages, in particular for blue supergiants, red supergiants, and for the immediate 
presupernova stage. 

We discuss the changes of the surface abundances due to rotationally induced mixing for main 
sequence and post main sequence stars. We single out two characteristics by which the effect 
of rotational mixing can be distinguished from that of massive close binary mass transfer, the 
only alternative process leading to non-standard chemical surface abundances in massive stars. 
A comparison with observed abundance anomalies in various types of massive stars supports the 
concept of rotational mixing in massive stars and indicates that it is responsible for most of the 
observed abundance anomalies. 

Subject headings: stars: rotation - stars: massive - chemical anomalies - surface properties 



1. Introduction 

Massive main sequence stars are rapid rota- 
tors. Equatorial rotation velocities span the range 
u ro t = 100 — 400 km s -1 , with B stars rotating clos- 
est to their break-up speed, w cr it (Fukuda 1982; 
Howarth et al. 1997). 

To date, stellar rotation and, in particular, ro- 
tationally induced mixing have not been taken into 
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account in extended studies of nucleosynthesis in 
massive stars. However, it is required to repro- 
duce, e.g., observed surface abundance anomalies 
of massive stars. In the present paper we predict 
the temporal evolution of the surface properties 
of rotating massive stars. This is of particular in- 
terest, as it allows a direct comparison of these 
models with observations. 

Besides abundance changes, lifetime and evo- 
lution in the Hertzsprung-Russell (HR) diagram 
of rotating stars differ from those of non-rotating 
stars. This is particularly interesting for age de- 
terminations of young stellar clusters. Since ap- 
propriate sets of rotating stellar models have not 
been available so far, these studies have mostly 
neglected effects of rotation. 

The calculations presented here treat rotation 



in stars according to the numerical method intro- 
duced by Endal & Sofia (1976, 1978); Pinsonneault 
et al. (1989). Our implementation is described in 
more detail in Heger et al. (2000). The results 
presented in this paper are based on calculations 
using the STERN code which contains a 36 iso- 
tope network including all stable isotopes up to 
silicon. As in Heger et al. (2000), we concentrate 
on Type II supernova progenitor stars, i.e., stars 
which evolve sufficiently massive cores to become 
core-collapse supernovae, but not massive enough 
to experience sufficient mass loss (or mixing) to be- 
come Wolf-Rayet stars (i.e., initial masses in the 
range ~ 8 — 20 M©). We also do not consider bi- 
naries. 

A list of our initial models and the evolution 
of stellar structure and angular momentum in the 
stellar interior are given in Heger et al. (2000). 
The discussion of the observable stellar properties, 
which is the topic of the present paper, is divided 
in two major parts: the main-sequence evolution 
(§2) and the post-main sequence evolution (§3). 
In §4 we compare our results to observations and 
give our conclusions in §5. 

2. Central hydrogen burning 

Core hydrogen burning comprises about 90 % of 
the total stellar lifetime and is therefore the evolu- 
tionary phase which can be most easily compared 
with observations. 

2.1. Nucleosynthesis and envelope abun- 
dances 

Central hydrogen burning in massive stars is 
governed by the CNO cycle (Bethe 1939; Clayton 
1968). Equilibrium mass fraction ratios of roughly 
1 : 50-100 : 1-5 for 12 C : 14 N : 16 O are established 
in the core (Arnould & Mowlavi 1993). This devi- 
ates distinctly from typical initial ratios of about 
1 : 0.3 : 3 (Grevesse & Noels 1993). Most of the 
initial carbon and oxygen is converted into 14 N. 
A subcycle of the CNO process, which transforms 
the initial 16 into 14 N, passes through 17 F and 
17 and is responsible for the production of 17 0. 
Further characteristic properties of the CNO cycle 
are the production of 13 C and the destruction of 
15 N. 

In the stellar layers above the convective core 
the CNO equilibrium abundances are not reached. 



In particular, there is almost no processing of hy- 
drogen into helium. However, isotopes which are 
produced at lower temperatures than where they 
are destroyed (e.g., 13 C and 17 0) can have abun- 
dance maxima in the envelope. 

Other nuclear cycles occurring during hydrogen 
burning involve, e.g., neon, sodium, magnesium, 
and aluminum. They do not contribute signifi- 
cantly to the energy generation, but they do al- 
ter the abundances of the involved isotopes. As 
the abundance changes of different isotopes occur 
at different temperatures, measuring their abun- 
dances at the stellar surface may contribute addi- 
tional information about the mixing processes in 
the stellar interior. 

The light isotopes 3 He, 6 Li, 7 Li, 9 Be, 10 B, and 
n B are destroyed at the low temperature in the 
stellar envelope. This allows to use them to con- 
strain mixing processes in the upper part of the 
envelope. For example, the surface depletion of 
lithium is applied to investigate models for the 
present sun (Talon & Charbonnel 1998). 

2.2. Internal mixing: the example of 
20 M© stars 

2.2.1. Mixing and barriers 

Fig. 1 compares the internal abundance profiles 
of the most abundant isotopes as a function of the 
Lagrangian mass coordinate, to, in a non-rotating 
and two rotating 20 M© models at core hydrogen 
exhaustion. The flat profiles in the innermost few 
solar masses are caused by convection. Small con- 
vective and/or semiconvective regions (Heger et al. 
2000) cause the "steps" in the profile above the 
convective core. 

Fig. 1 also shows the absence of mixing above 
this region in the non-rotating case, visible as a 
kink in the hydrogen, helium, and oxygen profiles 
at about to = 10 M©. Between to = 10 M© and 
to = 13 M©, a sharp decrease in 12 C occurs, ac- 
companied by an increase in 14 N. This is due to 
partial CNO processing "in place", i.e., without 
any transport. Deeper in the core, the carbon 
mass fraction increases again and reaches its equi- 
librium value. Below to — 10 M© 14 N approaches 
its CNO equilibrium value, while oxygen is de- 
pleted. There is no thermonuclear depletion of 
oxygen above to = 10 M©. Above to « 14 M©, the 
CNO isotopes retain their initial mass fractions. 
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In contrast, the rotating models show mixing of 
thermonuclear processed matter up to the surface 
of the star. In Model E20 (Fig. la; no inhibiting 
effects of gradients in the mean molecular weight 
(/i-gradients) on the rotationally induced instabili- 
ties) an extended helium gradient reaches from the 
upper edge of the convective core up to the sur- 
face. Along with helium, also nitrogen is enriched 
while carbon and oxygen are depleted. Between 
m = 5M Q and m = 10 M©, carbon is enhanced. 

If /i-gradients are allowed to inhibit the ro- 
tationally induced instabilities (Model E20B; 
Fig. 16), the /i-gradient which forms at the up- 
per edge of the convective core is not smoothed 
out, but instead almost completely chokes off any 
mixing between core and envelope quite early dur- 
ing the evolution. Above the "barrier" due to the 
/i-gradient (/i-barrier), however, mixing is efficient 
(Fig. 16). It's stronger than for Model E20, since 
the efficiency for compositional mixing was as- 
sumed to be f c = 1/30 in Model E20B instead 
of f c = 1/100 for Model E20 (Heger et al. 2000). 
Therefore, a stronger enrichment of nitrogen close 
to the surface results in Model E20B, due mainly 
to the processing of carbon into nitrogen. The 
16 O abundance in the envelope hardly changes 
in Model E20B, because the temperature above 
the /i-barrier is never high enough to allow for its 
processing. Helium, being the main reason for the 
/i-barrier, remains small in the envelope. 

In Fig. 2 we show the total diffusion coefficient 
as resulting from all mixing processes taken into 
account at the beginning, about the middle, and 
the end of central hydrogen burning for Models 
E15 and E15B. For Model E15B the /i-barrier can 
be seen a deep drop in the diffusion coefficient 
above the convective core at the two later times, 
while the /i-barrier is not yet present at the be- 
ginning of core hydrogen burning (solid line). In 
Model E15 such a /z-barrier never forms, since the 
rotationally induced mixing is assumed insensitive 
to /i- gradients, and thus material from the core can 
be mixed up to the surface during whole central 
hydrogen burning. However, note that the dif- 
fusion coefficient by rotationally induced mixing 
processes, dominated by Eddington-Sweet circu- 
lation in the envelope, is lower by a factor three 
compared to Model E15B, since f c is lower by this 
factor. Therefore, also the time-scale for mixing 
in the envelope is longer by that factor. 



2.2.2. Light element abundances 

The light elements, which are already destroyed 
at the rather low temperatures in the upper part 
of the envelope, allow to study the mixing close 
to the surface. Lithium is the most fragile and 
destroyed first, closely followed by beryllium. At 
somewhat higher temperatures, boron follows, and 
further in, more than five solar masses below the 
surface in the 20 M© models, the destruction of 
3 He proceeds (Fig. 3). 

In non-rotating stars, the isotopes are just de- 
stroyed "in place" wherever the temperature gets 
high enough, while at lower temperatures, i.e., 
closer to the surface, they are not altered (Fig. 3). 
Thus their surface abundance during central hy- 
drogen burning changes only if mass loss is suffi- 
ciently strong to remove the unaltered layers. This 
is the case in Model D20 for lithium and beryllium. 
Due to the strong temperature dependence of the 
destruction rates, the resulting abundance changes 
can become quite dramatic. Though Model D25 
experiences more mass loss than Model D20, its 
envelope is more extended and a more massive 
surface layer can preserve the initial stellar abun- 
dances till the end of central hydrogen burning - 
enough to result in lower depletion of lithium and 
beryllium then Model D20. 

In rotating stars, however, rotationally induced 
mixing causes a downward transport of the light 
elements. This has two important consequences. 
First, the light elements are transported into re- 
gions where the temperature is sufficiently high to 
destroy them. This leads to a depletion of these 
isotopes by the rotating massive stars during cen- 
tral hydrogen burning. The faster the mixing in 
the envelope, the stronger the depletion. Secondly, 
due to mass loss and expansion of the stars during 
central hydrogen burning, envelope layers which 
have been hot enough to destroy the light elements 
cool down, and the rotationally induced mixing 
can enrich these layers again with the residual of 
the light elements in the layers above before these 
outer layers get lost by winds. That is, some of 
the light isotopes can actually be preserved - to a 
small extent - by the rotationally induced mixing 
(e.g., Fig. 3a for lithium). We denote this effect 
as rotationally reduced depletion. 

Obviously, these two effects compete with each 
other. As long as mass loss is not large enough to 
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remove most of the unprocessed matter from the 
surface, both the surface abundance and the total 
mass of the corresponding isotope will be higher 
in the non-rotating star. But if mass loss is great 
enough to remove all the unprocessed matter of a 
non-rotating star, the rotating model can, for the 
same amount of mass loss, actually retain both a 
higher total mass and surface mass fraction of the 
fragile isotope. 

In the models of series "B" , where the mixing 
in the envelope is relatively fast, the destruction 
due to the mixing of the light elements into hot 
layers dominates. This effect increases with initial 
angular momentum, and it can destroy lithium, 
beryllium, and boron almost completely, and re- 
duce 3 He by more than a factor ten, even at the 
surface. 

2.2.3. Sodium and beyond 

Nuclear processing of the stable isotopes of 
neon, sodium, magnesium, and aluminum occurs 
only in the hot convective core of the star, i.e., be- 
low the /z-barrier. In the core, 21 Ne and 22 Ne are 
converted to 23 Na by the reaction 21 Ne (p,7) 22 Na 
(/3+^ c ) 22 Ne (p, 7 ) 23 Na. Thus a depletion of 21 Ne 
and 22 Ne is accompanied by an enrichment of 
23 Na. The long-lived radioactive isotope 26 Al is 
produced from 25 Mg by proton capture, but de- 
cays by electron capture to 26 Mg with a lifetime 
of ~ 1 Myr. The processing of 26 Mg to 27 Al leads 
to a slight enrichment of this aluminum isotope, 
but it is unimportant for 26 Mg. The abundances 
of 20 Ne and 24 Mg are essentially unaffected. 

No sodium or 26 Al enrichment was found in 
the envelope of the non-rotating Model D20, but 
strong abundance gradients of these species exist 
inside the envelope of the rotating Model E20 (in- 
efficient /i-barriers; see Fig. 5), and a slight enrich- 
ment is even found at the surface. In Model E20B 
the /x-gradient is weak enough to allow efficient 
transport of these processed isotopes from the core 
into the envelope only during the early phase of 
hydrogen burning. In particular, the rapid conver- 
sion of 22 Ne to 23 Na at hydrogen ignition allows it 
to be mixed into the envelope before the formation 
of the ^i-barricr. Due to the more efficient mixing 
inside the envelope, 23 Na is even more enriched 
at the surface than in Model E20, even though its 
total mass inside the envelope is much less. Sim- 
ilar processes occur for the other isotopes of the 



neon/sodium and the magnesium/aluminum cy- 
cles. 

The strength of enrichment of nuclear processed 
matter increases with initial mass (Heger et al. 
2000). Several effects contribute: With increasing 
stellar mass the stellar interior becomes increas- 
ingly dominated by radiation pressure, and the 
radiative zone of the star is closer to the adia- 
batic stratification, which, for example, increases 
the velocity of the Eddington-Sweet circulation 
(Heger et al. 2000). Additionally, the convective 
core comprises a larger fraction of the total stellar 
mass and mass loss increases. Both help to bring 
processed matter to the surface of the star. 

2.3. Evolution in the HR diagram 

At core hydrogen ignition, a rotating star be- 
haves, in some respect, as if it had less mass. It 
is less luminous (Heger et al. 2000), i.e., burns hy- 
drogen at a lower rate, and ages slower. This is 
due to the reduction of the effective gravity by 
the centrifugal force (Kippenhahn et al. 1970). At 
the same time, the radius of a rotating star is 
larger than that of a non-rotating star of same 
mass. Therefore, it appears in the HR diagram at 
lower effective temperature. This can be seen in 
Figs. 6, 7, and 8 by comparing the starting points 
of the evolutionary track of stars with different 
initial rotation rates in the HR diagram (Bodcn- 
hcimcr 1971; Endal & Sofia 1976). 

Without any rotationally induced mixing, ro- 
tating stars evolve essentially parallel to non- 
rotating stars in the HR diagram during core hy- 
drogen burning, but at lower luminosity and sur- 
face temperature (Kippenhahn et al. 1970; Endal 
& Sofia 1976). However, rotationally induced mix- 
ing increases the average mean molecular weight in 
the star by enriching the envelope in helium (§2.1; 
Heger et al. 2000). The degree of enrichment de- 
pends upon the inhibiting effect of ^-gradients and 
on initial angular momentum of the star (Heger 
ct al. 2000). In chemically homogeneous stars, 
L cx /j, 4 M 3 (Kippenhahn & Weigert 1991). That 
is, an increasing helium enrichment makes rotat- 
ing stars to evolve at higher luminosity which can 
exceed the luminosity of a non-rotating star of 
same initial mass (Fig. 6). The higher luminosity 
also causes the convective core to recede less or 
even to grow (Heger et al. 2000). 
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For faster initial rotation, the stars tend to live 
longer. At the beginning of hydrogen burning 
they are less luminous due to the reduced effective 
gravity, and use their fuel supply more sparingly 
(Fig. 10el-e3). Only at later times they evolve to 
higher luminosities because of the higher average 
mean molecular weight (Heger et al. 2000), but at 
the same time they have a larger fuel supply and 
larger convective cores. This latter effect is analo- 
gous to the effect of "convective core overshooting" 
(e.g., Schaller et al. 1992). The helium enrichment 
in the envelope can also cause the star to evolve 
to higher surface temperatures, and therefore to 
a position in the HR diagram where less evolved 
non-rotating stars of higher initial mass would be 
located (Langer 1992). 

These trends are demonstrated in Fig. 6 using 
the example of 12 M stars with different zero- 
age main sequence (ZAMS) rotational velocities 
and for the two cases, rotationally induced mix- 
ing being insensitive (Fig. 6a: Models D12, G12, 
E12, and F12) or sensitive (Fig. 66: Moddels D12, 
G12B, E12B, F12B, and H12B) to ^-gradients. In 
order to meet the observational constraints on the 
surface enrichment, the mixing efficiency is smaller 
in the former models (see Heger et al. 2000). 
In these models, the mixing occurs continuously 
without any hindrance by a /Lt-barrier (§2.2.1) and 
the stars gradually evolve towards higher luminosi- 
ties compared to non-rotating stars. They reach 
the end of central hydrogen burning at higher lu- 
minosities for higher ZAMS rotational velocity. 

As the fastest rotators remain close to uniform 
rotation, Model F12 hits the £!-limit of critical 
rotation (Langer 1997) at the end of central hy- 
drogen burning. Therefore, it can not contract 
and move back to higher temperatures before the 
"turn-off", like the other models shown in Fig. 6a. 
The same occurs for Model H15B. 

In the case where rotationally induced mixing 
is sensitive to ^-gradients, the envelope helium 
enrichment is only weak for the slowly rotating 
models (v Iot < 200 km s -1 ), since the /x-barrier be- 
comes effective early on. These stars evolve very 
similar to the non-rotating models. In the faster 
rotating models (see Models F12B and H12B in 
Fig. 66) the rotationally induced mixing initially 
operates efficiently enough to penetrate the grow- 
ing ^-barrier. The stars evolve to higher luminosi- 
ties and surface temperatures for some time and 



follow, in the HD diagram, a track almost perpen- 
dicular to the evolution of the non-rotating Model 
D12. During this time, their convective core mass 
remains about constant or even grows (in Model 
H12B). Eventually, the ^-barrier becomes strong 
enough to suppress the rotationally induced mix- 
ing, and the stars continue their evolution in the 
HR diagram parallel to that of the non-rotating 
model, but at a higher luminosity. This "^i-turn" 
occurs the later the higher the initial rotational 
velocity is (for a given initial mass). It also affects 
the evolution of the surface abundances (§2.5), be- 
cause it marks the point where the chemical evo- 
lution of the envelope becomes disconnected from 
that of the core. 

Very similar behavior is also found for the 
15 M (Fig. 7) and 20 M Q (Fig. 8) stars. The 
effects described above are stronger for the fast- 
rotating stars of lower mass because those are 
closer to their critical rotation rate. Additionally, 
the more massive stars are more dominated by ra- 
diation pressure in their interiors, and so helium 
enrichment affects them less. 

As a result of these dependences of the stel- 
lar evolutionary tracks in the HR diagram on the 
initial rotation rate, a given point in the HR di- 
agram is not uniquely related to a single initial 
mass, even for core hydrogen burning stars. In 
contrast, stars which differ in initial mass and ro- 
tation rate can evolve to the same position (Fig. 2 
of Flicgner et al. 1996). For example, the fast ro- 
tating Models H12B and H15B share, at the end 
of their central hydrogen burning, the same po- 
sition in the HR diagram as the slowly rotating 
models of 15 M Q and 20 M© (Figs. 66, 7, and 8). 
The rotating models can also reach positions in 
the HR diagram to the right of the main sequence 
band of non-rotating stars, i.e., the main sequence 
is broadened by the action of rotation. 

A distinguishing feature, however, between 
lower-mass fast rotators and higher-mass slow ro- 
tators which share the same position in the HR 
diagram, is the different evolutionary time they 
need to get to that position. This can be seen 
in Fig. 10el-e3: For example, the fast rotating 
Model H12B reaches the same luminosity and also 
about the same surface temperature after 25 Myr, 
while the slowly rotating 15 M© stars terminate 
central hydrogen burning at the same position 
already after only 10 Myr. Another such distin- 
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guishing feature is the expected surface abundance 
pattern, which is discussed in §2.5. 

Hot stars may become Wolf-Rayet stars with 
surface hydrogen mass fractions below 40%, ei- 
ther due to mass loss or due to mixing (Maeder 
1982; Langer 1987). As mentioned above, stars 
with more internal mixing also have more mass 
loss due to their higher luminosity. As an exam- 
ple, the extremely rapidly rotating Model H20B 
becomes a Wolf-Rayet star at a central hydro- 
gen abundance of 4%, i.e., during central hydro- 
gen burning (Fig. 8). Thus, rotationally induced 
mixing can significantly influence the lower limit- 
ing initial mass required for Wolf-Rayet star for- 
mation during central hydrogen burning in single 
stars (Maeder 1987, 1999; Langer 1992; Fliegner 
& Langer 1995). 

2.4. Rotational velocities 

The evolution of the surface rotational velocity 
for various models is displayed in Fig. 10dl-d3. It 
is affected by several different contributions: First, 
the loss of angular momentum due to stellar winds 
leads to a slow-down of the star. For the models 
shown in Fig. 10dl-d3, this reduces the initial an- 
gular momentum of the 20 M models by < 40 %, 
by < 20% for the 15 M models, and < 10% for 
the 12 M models. The second effect is the in- 
crease of the total moment of inertia of the star, 
for the stars considered here by about a factor two 
from core hydrogen ignition until shortly before 
core hydrogen exhaustion, due the expansion of 
the envelope. Both the moment of inertia and the 
angular momentum of the models are dominated 
by the envelope, which stays close to uniform ro- 
tation. Third, the increase of the stellar radius by 
about a factor two to three during core hydrogen 
burning leads to an increase of the surface rota- 
tional velocity of the stars by the same factor for 
a given angular velocity, v mt = rui. This effect al- 
most completely cancels out the slowing down due 
to the two preceding effects, as can be seen from 
the small change in the surface rotational velocity 
in Fig. I0dl-d3. 

The rise in rotational velocity at the end of core 
hydrogen burning is due to the overall contraction 
of the star before ignition of hydrogen shell burn- 
ing. The deep drop of the rotational velocity in 
Model H20B (Fig. 10<£2) is caused by the angu- 
lar momentum loss occurring as a consequence of 



the strong mass loss as it becomes a Wolf-Rayet 
star. Model H12B (Fig. lOdl ) reaches the O-limit 
(Langer 1997) at the end of core hydrogen burn- 
ing, and therefore strong anisotropic mass loss is 
expected in this case (Langer 1998). 

2.5. Surface abundances 

The efficiencies of the rotationally induced mix- 
ing in our models have been calibrated (Heger 
et al. 2000) such that the average surface abun- 
dance trends for helium and CNO abundances on 
the main sequence are consistent with the obser- 
vations. However, the enrichment or depletion of 
these elements as a function of rotation rate, time, 
and initial mass were not adjusted. Also, the sur- 
face abundances of the rare helium and CNO iso- 
topes, as well as the abundances of all other ele- 
ments were not considered in the calibration pro- 
cess. It is thus essential for a comparison with 
observations to also investigate the surface abun- 
dances of our models during the main sequence 
phase. 

2.5.1. 3 He, Li, Be, B 

An important feature of rotationally induced 
mixing is the depletion of light elements, like 
lithium and boron, at the surface. These ele- 
ments are destroyed at much lower temperatures 
than those necessary for CNO processing. There- 
fore, much stronger abundance changes are in- 
duced by rotation. The rotating main sequence 
models computed by Fliegner et al. (1996) pre- 
dicted that stars which show enrichment in nitro- 
gen should be strongly depleted in boron, which 
was found in good agreement with the observa- 
tional data by Venn et al. (1996). This finding is 
supported by the present work (Figs. 9al-a3 and 
llal-a3): the boron depletion is much stronger 
and appears earlier than the nitrogen enrichment. 
Even the slow rotators with ZAMS rotational ve- 
locities of ~ 100 km s -1 , which barely show any 
nitrogen enrichment, become depleted in boron by 
- 0.5 dex. 

At the end of central hydrogen burning, the 
non-rotating Model D20 shows a lithium deple- 
tion which is more than two orders of magnitude 
greater than that of the slowly rotating Model 
G20. This is due to the rotationally reduced deple- 
tion (§2.2.2). In a similar way, Model E20, where 
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the efficiency of the mixing in the envelope is less 
than in Model E20B, experiences less destruction 
of lithium than Model D20 (Fig. 3a). But the de- 
struction is also less than in Model E20B, where 
the efficient mixing would lead to a strong deple- 
tion even without mass loss (Fig. 36; see also Fig. 4 
for a model with less mass loss). 

While beryllium is preserved at the surface in 
most of our non-rotating models, except for Model 
D20, already our slowest rotating models, Models 
G12B, G15B, and G20B, show considerable de- 
struction of beryllium at the end of central hy- 
drogen burning. Similar to lithium, those models 
where the inhibiting effect of the ^-gradients on ro- 
tationally induced mixing are taken into account 
but a more efficient chemical mixing is assumed, 
show much stronger depletion. Due to the stronger 
mass loss, the depletion increases with increasing 
initial stellar mass, exceptions being, again, Model 
D20, and Models H12B, H15B, and H20B which 
show an opposite trend, but at very high depletion 
factors (Table 1). 

The 3 He isotope is much less fragile than the 
species discussed above and is not depleted at the 
surface in any of the non-rotating models and only 
slightly depleted in the slow to average speed rota- 
tors (< 200 kms -1 ). In the fast rotators, we find 
a depletion of up to almost a factor 100 (Model 
H20B; Table 1). The magnitude of the depletion 
depends more strongly on the initial rotation than 
on the initial mass. 

2.5.2. C, N, O 

The evolution of the surface abundances of 12 C, 
14 N, and 16 relative to their initial values, and 
that of the 4 He mass fraction during central hydro- 
gen burning are shown in Fig. 9 for 12 M Q , 15 M Q 
and 20 M Q stars with different initial angular mo- 
menta. Fig. 10 displays the corresponding ratios of 
elemental carbon, nitrogen and oxygen abundance 
relative to their initial values. The surface abun- 
dances at core hydrogen exhaustion are shown in 
Table 1 for elements up to sodium. 

The dominant isotope in the CNO cycle, 14 N, 
is enriched at the surface at the expense of 12 C 
and 16 0. This enrichment is very small for the 
models with a ZAMS equatorial rotational velocity 
of <~ 100 km s^ 1 (< 0.1 dex), but for models with 
typical rotational velocity (~ 200 kms -1 ) the ni- 



trogen enrichment at the end of central hydrogen 
burning is a factor two to three, and it is increasing 
with initial mass. At higher initial rotation rates 
the final enrichment factor is about 6 and 10 for 
the models with ~ 300 km s -1 and ~ 400 km s -1 
ZAMS rotational velocity, respectively. The latter 
value is already close the CNO equilibrium, i.e., 
most of the initial 12 C and 16 are converted into 
14 N. A notable depletion of 12 C occurs only for 
ZAMS rotational velocities of > 200 kms -1 , but 
becomes as high as a factor ~ 30 for the fastest ro- 
tating models — which is even higher than the de- 
pletion factor of 12 C in the stellar center, since the 
12 C CNO equilibrium abundance is lower for lower 
temperatures (Arnould & Mowlavi 1993). The de- 
pletion of 16 and the enrichment of 4 He are cor- 
related and become significant only for stars with 
ZAMS rotational velocities of > 300 kms -1 , the 
stronger the higher the initial mass of the stars. 
The coupling of these two isotopes is due to the 
/x-barrier: both, the destruction of 16 and the 
production of 4 He occur within or below the /x- 
barrier (Heger et al. 2000). 

In the extremely rapidly spinning 20 M Model 
H20B the surface helium enrichment becomes so 
large (> 60% by mass), at a central hydrogen 
mass fraction of ~ 4 %, that the star can be consid- 
ered a Wolf-Rayet (WR) star, with a correspond- 
ingly high mass loss rate. Even though the star 
experiences high mass loss during that WR phase, 
it becomes a red supergiant after core hydrogen 
exhaustion. The high mass loss rate reduces the 
rotational velocity and the luminosity of the star, 
as shown in Figs. 10d3 and 10e5. 

The strength of the surface enrichment is cor- 
related with the time when notable enrichment 
first appears during the evolution. This is a sign 
of an over-all increasing strength of the mixing 
with increasing initial angular momentum. It is 
interesting to note that the /x-turn in the HR di- 
agram (§2.3) for the fast rotating models also co- 
incides to a "kink" in the evolution of the surface 
abundances, especially for isotopes produced be- 
low the /x-barrier, namely 16 and 4 He, and, to 
some extent, 14 N. These /i-turns appear at an 
age of about lOMyr, 20Myr, 15Myr, and lOMyr 
for Models F12B, H12B, H15B, and H20B, respec- 
tively (Fig. 9). 
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2.5.3. F, Na, Al, Mg 

Other interesting changes in the surface abun- 
dances due to rotationally induced mixing are 
the depletion of fluorine ( 19 F) and enrichment in 
sodium ( 23 Na), both the only stable isotopes of 
these elements (Figs. Ilbl-b3,cl-c3). An impor- 
tant product of the magnesium-aluminum cycle 
(Arnould & Mowlavi 1993) is the production of 
the radioactive isotope 26 Al by 25 Mg (p,7) 26 Al in 
the central convective core. Typically it reaches a 
mass fraction of a few 1CP 6 in the center for the 
initial composition used in the present work. This 
isotope has a half life of 0.716 Myr and emits 7- 
rays due to its /3-decay, which has been observed 
in the interstellar medium of our Galaxy (Ober- 
lack et al. 1996). The surface abundance of 26 Al 
shows the interplay of production, decay and mix- 
ing (Figs. Ildl-d3): a steep initial increase in the 
abundance appears when 26 Al is mixed to the sur- 
face within the first few Myr. For the fast rotat- 
ing models, the abundance slowly continues to in- 
crease until the mixing between core and envelope 
is shut off at the /x-turn, and the further evolu- 
tion of the 26 Al surface abundance is mainly gov- 
erned by its decay. The drop of the 26 Al surface 
abundance, as, e.g., seen in Fig. lldl , corresponds 
to an exponential decay. The apparently different 
slopes in Figs, \ldl-d3 are due to the different 
time axes. 

The surface abundance of 26 Al strongly in- 
creases with the amount of mixing, i.e., with the 
ZAMS rotational velocity of the star. But it is 
also a sensitive function of the initial mass: The 
fastest rotating 20 M Q star (Model H20B) has 
an about two orders of magnitude higher surface 
abundance than the fastest rotating 12 M star 
(Model H12B). This difference is even larger for 
the slower rotators. 

The larger surface enrichment in 26 Al in the 
more massive stars is due mostly to their shorter 
lifetime. More massive stars show a comparable 
enrichment in products of central hydrogen burn- 
ing as a function of the fraction of their lifetime, 
within a factor two of less massive stars, as can be 
seen, e.g., from the surface enrichment in helium 
(Fig. 9dl-d3). But during the shorter absolute life- 
time, less 26 Al decays. 

Several interesting changes in the surface iso- 
topic ratios occur for different elements as well. 



They are shown in Heger (1998). In contrast to 
the changes in the elemental abundances, these 
isotopic ratios are not easily measurable in the at- 
mospheres of hot stars, and therefore cannot be 
used as observational tests. 

3. Post-central hydrogen burning 
3.1. Evolution in the HR diagram 

After core hydrogen exhaustion the stars dis- 
cussed in the present work become red super- 
giants. They evolve to the the Hayashi line within 
a few lOOOOyr. In Figs. 12 and 13, the evolu- 
tionary tracks in the HR diagram of the rotat- 
ing models with a ZAMS rotational velocity of 
<~ 200 km s -1 and non-rotating models are com- 
pared. The tracks go from the ZAMS up to core 
collapse, except for the lowest mass Models D10, 
D12, and E08 in Fig. 12, which end before neon ig- 
nition. The evolution of Model E25 is only shown 
until it becomes a Wolf-Rayct star during central 
helium burning, where it loses its hydrogen-rich 
envelope. In Figs. 14, 15, and 16 the evolution of 
stars with different initial rotation rates are com- 
pared. 

Since the rotating stars develop larger helium 
cores, they become more luminous red supergiants 
than the non-rotating stars of equal mass. This 
is shown in Fig. 12 for models with inefficient /i- 
gradients, and in Fig. 14, for different initial ro- 
tational velocities. In the former case, the models 
also evolve to notably lower effective temperatures 
than their non-rotating counterparts. These mod- 
els have considerably larger helium cores and, due 
to their stronger mass loss, a by far higher enve- 
lope to core mass and mass to luminosity ratio. 
Due to this, these stars, especially the most lu- 
minous of them, become highly unstable to radial 
pulsations (Heger et al. 1997). 

In models where the effect of /u-gradients on ro- 
tationally induced mixing is taken into account, 
the effective temperatures reached on the Hayashi 
line are not significantly different from those of the 
non-rotating models. The increase in luminosity is 
much less (Figs. 15 and 16), except for the fastest 
rotators, like Model H12B (see Heger 1998). This 
model does not become a red supergiant immedi- 
ately after core hydrogen exhaustion, but burns 
most of the helium as a blue supergiant before it 
moves to the Hayashi line. 
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Stars with initial masses of < 12 Mq show so- 
called "blue loops", i.e., they evolve from red into 
blue supergiants (Figs. 12, 13 and 14). Stars with 
more rotationally induced mixing evolve to lower 
effective temperatures and to higher luminosities 
during their blue loops (Figs. 12, 13, and 14). The 
highly peculiar evolution of the surface rotation 
rate and mass loss during a blue loop is discussed 
in (Hcgcr & Langer 1998). Here we give in Ta- 
ble 4 the average radii rotational velocities and 
specific angular momenta during the blue loop. 
The maximum rotation rate, which usually is ob- 
tained briefly after the transition to a blue su- 
pergiant, can be much higher than that average 
(§3-4). 

3.2. First dredge-up abundances 

Red supcrgiant have an extended surface con- 
vection zone which comprises the major part of the 
hydrogen-rich envelope. Since convective mixing is 
much faster than rotationally induced mixing, the 
convective envelope is essentially chemically ho- 
mogeneous. The thermonuclear processed matter 
which becomes engulfed into the convective en- 
velope appears "instantaneously" at the surface, 
compared to the evolutionary time-scale of the 
star. This rapid mixing of processed matter to the 
surface, which occurs also in non-rotating stars, is 
known as a dredge-up, and its first occurrence af- 
ter core hydrogen depletion is accordingly called 
the first dredge-up. 

In Fig. 17 thin lines show, as a function of the 
ZAMS rotational velocity, the surface mass frac- 
tion of 12 C, 14 N, and 16 relative to their initial 
values and that of 4 He at core hydrogen exhaus- 
tion for three different initial masses for models 
where rotationally induced mixing was assumed 
to be sensitive to /i-gradients. Thick lines give 
the same quantities after the first dredge-up. The 
surface abundances of elements up to sodium after 
the first dredge- up are also given in Table 2. 

The slowly and non-rotating models show a sig- 
nificant enhancement in 14 N due to the dredge-up. 
In contrast, in the fast rotators the changes of the 
surface abundances due to dredge-up are small in 
comparison to the surface enrichment which al- 
ready occurred during central hydrogen burning. 
That is, the envelope is already well-mixed before 
the dredge-up. 



Since, during the dredge-up, the convective en- 
velope penetrates into regions which were located 
below the ^-barrier, the surface carbon mass frac- 
tion is reduced by diluting the 12 C in the envelope 
with the 12 C-dcpleted matter from deeper inside 
the star. In a similar way, 16 is depleted and 4 He 
is enriched. The bottom of the convective enve- 
lope does not get hot enough for nuclear burning, 
i.e., the abundance changes at the surface during 
the dredge-up result only from mixing of already 
processed matter into the convective envelope as 
it extends to lower mass coordinates. 

In the case of rotationally induced mixing, con- 
version of 12 C into 14 N occurs above the /z-barrier 
and therefore the composition of matter originally 
located close to the surface of the stars can be al- 
tered by first mixing it down, then processing it, 
and mixing it up again. Due to this, the surface 
enrichment in the faster rotating stars is stronger, 
both before and after the dredge-up, than in the 
slowly or non-rotating stars. The stronger deple- 
tion of 16 and enrichment of 4 Hc for the faster 
rotating stars is mainly because in those the ro- 
tationally induced mixing can overcome the fi- 
barrier for a longer time during core hydrogen 
burning than in the slow rotators. 

3.3. Second dredge-up abundances 

After core helium exhaustion, the bottom of the 
convective envelope of the red supergiants moves 
to lower mass coordinates. Thus, it can trans- 
port matter which has been processed by the hy- 
drogen shell to the surface of the star. However, 
in most of the models considered here, the con- 
vective envelope does not penetrate much deeper 
during the second dredge-up than it did during 
the first dredge-up. Additionally, not much mix- 
ing occurs in the radiative part of the envelope 
during core helium burning, and the temperatures 
are mostly too low for significant thermonuclear 
processing. Thus, with the exception of the low- 
est masses (e.g., Models D10 and E08), no strong 
changes in the surface abundances occur during 
the second dredge-up. In stars which experience a 
blue loop, the processing can be slightly stronger, 
since during this loop the envelope contracts and 
a larger part of it can become sufficiently hot for 
thermonuclear reactions to occur. 

The surface abundances of the most important 
isotopes after the second dredgc-up are given in 
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Table 3. Since in the case of the second dredge- 
up the mixing in the convective envelope is also 
fast enough to make it chemically homogeneous, 
Table 3 gives the presupernova abundances of the 
entire envelope. 

3.4. Rotational velocities 

The majority of the stellar models investigated 
here expand to red supergiants after core hydrogen 
depletion. The expansion of the stellar envelope 
increases the moment of inertia and the angular 
velocity of the rigidly rotating convective envelope 
of the star decreases considerably. This is the most 
important effect contributing to the slow-down of 
the star's rotation when evolving from the main 
sequence to the red supergiant regime. Addition- 
ally strong mass loss sets in during the red super- 
giant phase, and the corresponding loss of angular 
momentum due to stellar winds becomes impor- 
tant. This is also true for the models of lower 
initial mass, in which angular momentum loss was 
almost negligible during central hydrogen burn- 
ing. Since the rise of the moment of inertia is 
much larger than the growth in radius (Heger & 
Langer 1998), the surface rotational velocity de- 
creases considerably and drops below lkms -1 in 
most models (Table 4, Fig. 18). During the second 
dredge- up, i.e., after core helium exhaustion, the 
models become even more luminous and extended 
(Table 4). 

Models with an initial mass of 12 M Q or be- 
low, experiencing a blue loop, show a remarkable 
spin-up during this phase and can lose the major 
part of their total angular momentum (see Heger 
& Langer 1998). This explains the small difference 
in the angular velocity of the 12 M Q models com- 
pared to the corresponding 15 M Q models, which 
do not develop blue loops, in relation to the dif- 
ferences between other models with same input 
physics but different initial mass (see e.g., Models 
E12B, E15B, and E20B in Table 4). 

Blue supergiants born out of such a transition 
can rotate close to their critical rotation rate. This 
fast rotation occurs because the envelope angular 
momentum is concentrated in a small mass frac- 
tion at the surface of the star. Due to mass loss the 
rapidly rotating layers get lost soon and the star 
slows down quickly (see Heger & Langer 1998). 

Model H12B evolves as a blue supergiant during 



the first part of core helium burning. Since it has 
a large initial angular momentum and does not ex- 
perience any significant previous angular momen- 
tum loss, it evolves close to critical rotation during 
this phase. Even the enhancement of the mass loss 
rate due to rotation does not cause a sufficiently 
high angular momentum loss to slow it down from 
critical rotation. 

Models E25 and H20B become Wolf-Rayet stars 
during central helium burning. The evolution of 
Model E25 has been followed until core collapse. 
The strong mass loss during the Wolf-Rayet phase 
rapidly slows down the rotation of the helium en- 
velope. None of the instabilities considered in the 
present work is fast enough to keep it in rigid ro- 
tation with the core or the rest of the helium-rich 
layers. Thus strong differential rotation occurs in 
the envelope and the core is not slowed down sig- 
nificantly. The slower rotation of the core com- 
pared to other models (see Heger et al. 2000) is 
mostly due to the loss of angular momentum dur- 
ing core hydrogen burning. 

4. Observational evidence 

There is a considerable body of observations of 
abundance anomalies in massive stars. In par- 
ticular, helium and CN-abundances of main se- 
quence O and B stars by Gies & Lambert (1992); 
Herrero (1994); Vrancken et al. (2000) have been 
used to calibrate the mixing efficiencies for our 
models (Heger et al. 2000), which thus reproduced 
this on average. For a moderate initial rotation 
rate, our models obtain a rather insignificant he- 
lium enrichment (~1% by mass), a carbon deple- 
tion of less than 50 %, and a nitrogen enhancement 
of more than a factor two (Table 1). Although the 
observational facts about OB main sequence stars 
are not completely unambigous, the quoted num- 
bers appear to be supported by recent empirical 
studies (McErlean et al. 1998; Smith & Howarth 
1998; Herrero et al. 1999; Korn et al. 1999; McEr- 
lean et al. 1999). However, our main sequence 
models provide the following testable predictions: 

1. Stronger surface anomalies should occur for 
larger rotation velocities, v Iot . A test of this 
feature is complicated by two things, namely 
that not v rot but only v rot sin i is accessible 
to observations, where i is the inclination an- 
gle, and that abundance determinations are 
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often only possible for stars with relatively 
small values of u ro t sin i (cf. Gies & Lambert 
1992). Nevertheless, some tentative support 
for a positive correlation of enrichment and 
rotation comes from McErlcan ct al. (1999) 
and Herrero et al. (1999), although other in- 
terpretations of the observational data ap- 
pear possible as well. 

2. Stronger surface anomalies should occur for 
more luminous stars. This seems in fact to 
be the case, according to McErlean et al. 
(1999) and Gies & Lambert (1992). 

3. Lithium, beryllium and boron are consider- 
ably depleted. While there is no observa- 
tional data on the first two elements for 
hot luminous stars, the boron abundance in 
B star has been investigated by Venn et al. 
(1996), and by Cunha et al. (1997). Fliegner 
et al. (1996) interpreted the boron depletion 
in B stars which at the same time are not en- 
riched in nitrogen as evidence for rotational 
mixing. 

There is also evidence from more evolved lu- 
minous stars supporting the scenario of rotational 
mixing of the kind obtained by our models. How- 
ever, an additional problem appears in this case, 
namely that — for stars in the mass range con- 
sidered here — the post main sequence radius 
or temperature evolution can not yet be reliably 
predicted by stellar evolution theory (Langer & 
Maeder 1995). Therefore, it is often not unam- 
biguous which fraction of an abundance anomaly 
in an evolved massive star might be attributed 
to a possible previous convective dredge-up which 
would have occurred were the star a red supergiant 
in the past. 

However, Venn (1995a, b, 1999) has analyzed 
abundances in Galactic and Magellanic Cloud 
A supergiants and found results which are incom- 
patible with standard evolution with or without 
red supergiant dredge-up but rather call for an 
additional enrichment process like rotationally in- 
duced mixing. 

Further abundance anomalies which might be 
related to rotational mixing are the strong sig- 
nature of CNO processing in the progenitor of 
SN 1987A (Frannson ct al. 1989), small 12 C/ 13 C- 
ratios in red giants (Briley et al. 1995), sodium 



enrichment in cool supergiants (Boyarchuck et al. 
1988; Hill 1997; Hill et al. 1997), and simultaneous 
carbon and nitrogen enrichment in several Wolf- 
Rayet stars (Crowthcr et al. 1995; Langer 1991, 
for more details). 

While perhaps none of the quoted observational 
evidences is completely convincing by itself, their 
huge amount gives strong support for the idea that 
a non-standard (i.e., non-convective) enrichment 
process must in fact operate in massive stars. As 
the rotational mixing during core hydrogen burn- 
ing can also increase the luminosity-to-mass ra- 
tio, the so called mass discrepancy, i.e., the find- 
ing that masses of massive main sequence stars 
derived through spectroscopy or by constraints 
from radiation driven wind models are on aver- 
age smaller than masses derived from stellar evolu- 
tionary tracks (cf. Herrero et al. 1992, 1998), gives 
further support to extra mixing in massive stars. 

The only scenario presently competing with the 
rotational mixing scenario is that of close binary 
mass transfer. In fact, many massive stars are 
members of close binary systems, and one should 
expect mass transfer to contribute to the observed 
abundance anomalies in luminous stars. In a re- 
cent comprehensive study of massive close binary 
evolutionary models, Wcllstcin & Langer (1999) 
and Wellstein et al. (2000) find in fact that the 
enrichment of the accreting star in a binary may 
closely resemble the enrichment produced by ro- 
tational mixing in the models presented here. In 
particular, the CNO pattern, and the helium and 
sodium enhancements are found to be very similar. 
Moreover, those binary components which accrete 
significant amounts of mass from their companion 
end up in rather wide orbits, with companions of 
much lower mass. That is, after the supernova 
explosions of the companions, the accretion stars 
may not be easily recognizable by a large space 
velocity. Due to the accretion of angular momen- 
tum they might also be rapid rotators, which may 
lead to a positive correlation of enrichment and 
rotational velocity. 

Nevertheless, we identify two distinguishing 
characteristics of the rotation and the binary en- 
richment scenarios. First, as already anticipated 
by Fliegner et al. (1996), boron depleted main se- 
quence stars which are not nitrogen enriched are 
only produced by rotation. While Fliegner et al. 
(1996) used estimates to exclude such abundance 
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pattern for the binary scenario, this is now con- 
firmed in detail by the recent grid of binary models 
of Wellstein et al. (2000). Second, the rotational 
mixing scenario predicts a stronger enrichment for 
more luminous stars (see above). In contrast, the 
binary models of Wellstein & Langer (1999) and 
Wellstein et al. (2000), which cover a luminosity 
range of logL/L Q ~ 4.5 .. . 6.0, do not show this 
trend. 

As there is empirical evidence for both distin- 
guishing predictions of the rotational mixing sce- 
nario (cf. Fliegner et al. 1996; McErlean et al. 
1999), the contribution of the binary scenario 
seems not to be the dominant one. However, more 
observations are clearly needed to prove its signif- 
icance. We suggest that in particular the com- 
bined study of boron and CNO elements in hot 
stars will allow to draw more stringent conclu- 
sions, and that a simultaneous derivation of the 
rotational and also of the radial velocities would 
be very valuable to provide further clues as to the 
relative importance of the mass transfer and the 
rotational mixing scenario for an understanding of 
the multitude of abundance peculiarities in mas- 
sive stars. 

5. Conclusions 

Observations show that evolved massive main- 
sequence stars have nuclear processed matter at 
their surfaces (§4). Our models show that these 
surface abundance anomalies can be well explained 
by rotationally induced mixing. While the mix- 
ing processes we consider are gauged to repro- 
duce the average surface enrichment of nitrogen 
and helium, they also reproduce distinctive abun- 
dance trends that are not expected in competing 
binary scenarios, where the surface anomalies arc 
due to mass transfer from the companion. Specifi- 
cally, our models produce a boron depletion with- 
out nitrogen enrichment and an increasing nitro- 
gen abundance with increasing luminosity. 

Due to the strong temperature dependence of 
the destruction reactions, light isotopes will essen- 
tially be either preserved or completely destroyed 
at the surface of non-rotating stars, depending on 
the strength of the mass loss. Rotationally in- 
duced mixing, on the other hand, leads to a grad- 
ual depletion. Rotationally reduced depletion can 
also preserve some of the light elements in cases 



where a non-rotating stars would show complete 
depletion, although the resulting depletion might 
still be quite considerable (Table 1). 

Due to rotationally induced mixing, the chem- 
ical signature of hydrogen burning is apparent at 
the stellar surface already on the main sequence. 
Prominent features are enrichments of 4 He, 13 C, 
14 N, 17 , 23 Na, and 26 Al, while 12 C, 15 N, 16 0, 
18 0, and 19 F are depleted (Table 1). The faster 
the initial rotation the stronger the resulting mix- 
ing and the longer the mixing between core and en- 
velope can be maintained until a sufficiently large 
/i-barrier forms which suppresses further mixing 
(Fig. 11). 

The mixing between core and envelope in- 
creases the mean molecular weight of the star 
and leads to larger convective core and higher lu- 
minosity at the end of central hydrogen burning, 
i.e., the main-sequence band is broadened (Figs. 12 
and 13). The larger fuel supply, but even more so 
the lower initial luminosity of the rotating stars, 
helps to prolong the time they spend on the main 
sequence (Fig. 10). As the track in the HR di- 
agram can be significantly changed by rotation, 
a given position on the main sequence band does 
not unambiguously correspond to one initial mass 
(§2.3), i.e., there is no unique mass- luminosity re- 
lation for main sequence stars. Despite mass loss 
and expansion, the stellar rotation velocity stays 
about constant on the main sequence (Fig. 10). 

The first dredge-up changes the surface abun- 
dances by a factor which is smaller the larger the 
initial rotation rate is, since the fast rotators are 
already well mixed at the end of central hydro- 
gen burning (Fig. 17). Even after the dredge-up, 
the fast rotators maintain chemical peculiarities, 
as, e.g., a significantly stronger CNO processing. 
The second dredge-up does hardly affect the sur- 
face abundances, except in the two low-mass cases 
where the surface convection zone penetrates into 
the helium shell. Some small abundance changes 
are also present in the models that undergo a blue 
loop, since there the hydrogen shell burning is 
more extended in mass. However, after the blue 
loop, these red supergiants are rotating slower by a 
factor ~2-3, compared to red supergiants without 
a previous blue loop. 
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Table 1: Surface abundances at core hydrogen exhaustion for the indicated model sequences. Given are the 4 He mass fraction, and the 
logarithm of the ratios of the 3 Hc, lithium, boron, beryllium, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 0, 19 F, and 23 Na mass fractions relative to 
their initial values. 
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Table 2: Surface abundances after the first dredge-up for the indicated model sequences. Given are the 4 He mass fraction, and the logarithm 
of the ratios of the 3 He, lithium, boron, beryllium, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 0, 19 F, and 23 Na mass fractions relative to their initial 
values. 
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Table 3: Surface abundances after the second dredge-up for the indicated model sequences. Given are the 
helium mass fraction, and the logarithm of the ratios of the lithium, boron, beryllium, nitrogen, carbon, 
oxygen, fluorine, and sodium mass fractions relative to their initial values. 
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Table 4: Radius, R, equatorial rotation velocity, v, and equatorial specific angular momentum, ji% (in units of 10 18 cm 2 s^ 1 ) at the ZAMS, 
during the RSG phase (average values), during the blue loop (average values), and for the final models. Additionally, maximum values of 
rotation velocity, v max , and specific angular momentum, jis.max, during the blue loop are given. Dashes indicate the absence of a blue loop. 
Lifetime during the MS, RSG, and BSG phases are respectively given in columns tms, t rsg, and tbsg- 




Fig. 1. — Mass fractions of different isotopes as a function of the mass coordinate, m, at core hydrogen 
exhaustion. Compared are the chemical structures of rotating (thick lines) and a non-rotating (thin lines; 
same in both Panels; Model D20) 20 M Q models. The rotating models have a ZAMS equatorial rotational 
velocity of ~ 200 krns" 1 . Panel a: Model E20, where rotationally induced mixing is not inhibited by \x- 
gradients. Panel b: Model E20B, where rotationally induced mixing is (slightly; / M = 0.05) inhibited by 
/x-gradients. See also Figs. 3 and 5. 
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Fig. 2. — Total diffusion coefficient for chemical mixing as a function of mass coordinate, to, at three 
different stages of central hydrogen burning: 1 % of central hydrogen, by mass, burnt (X c = 0.69; solid 
line), and at central hydrogen mass fractions of 30% (X c = 0.30; dotted line) and 1% (X c — 0.01; dashed 
line). Panel a: Model E15, rotationally induced mixing is not inhibited by /x-gradients. Panel b: Model 
E15B, rotationally induced mixing is (slightly; — 0.05) inhibited by ^-gradients. The labels in the figure 
indicate the regime of diffusion coefficient for convection and Eddington-Sweet circulation (ES circ). For 
Model E15B, the /x-barrier is indecated. In that region, the Eddington-Sweet circulation is inhibited and 
the Goldreich-Schubert-Fricke instability (GSF) is most efficient process. Note the region of essentially no 
mixing above the central convection zone for the two evolved stages (X c = 0.30, 0.01) of Model E15B. 
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Fig. 3. — Mass fractions of 3 He (solid lines), lithium (dotted lines), 9 Be (dashed lines), boron (dash-dotted 
lines), and 19 F (dash-triple-dotted lines) as a function of the mass coordinate, m, at core hydrogen exhaustion 
for the same models as in Fig. 1. See also Fig. 5. 
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Fig. 4. — Mass fractions of 3 He (solid lines), lithium (dotted lines), 9 Be (dashed lines), boron (dash-dotted 
lines), and 19 F (dash-triplc-dottcd lines) as a function of the mass coordinate, m, at core hydrogen exhaustion. 
Compared are the chemical structures of rotating (thick lines) and a non- rotating (thin lines; same in both 
Panels; Model D12) 12 M Q models. The rotating models have a ZAMS equatorial rotational velocity of 
~ 200 kms -1 . Panel a: Model E12, where rotationally induced mixing is not inhibited by /i-gradicnts. 
Panel b: Model E12B, where rotationally induced mixing is (slightly: / M = 0.05) inhibited by ^-gradients. 
See also Fig. 3. 
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Fig. 5. — Mass fractions of 22 Ne (solid lines), 23 Na (dotted lines), 25 Mg (dashed lines), and 26 Mg (dash- 
dotted lines) as a function of the mass coordinate, m, at core hydrogen exhaustion for the same models as 
in Fig. 1. See also Fig. 3. 
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Fig. 6. — Main sequence evolution of 12 M Q stars with different ZAMS rotational velocities in the HR 
diagram. Panel a: Models D12, G12, E12, and F12, where rotationally induced mixing is assumed insensitive 
to /j-gradicnts. Panel b: Models D12, G12B, E12B, F12B, and H12B, where rotationally induced mixing is 
assumed to be sensitive to ^-gradients (/ M = 0.05). 
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Fig. 7. — Core hydrogen burning evolution in the HR diagram of 15 M Q models with ZAMS rotational 
velocities of Okm s" 1 (Model D15; solid line), 102kms- 1 (Model G15B; dotted line), 206 kms" 1 (Model 
E15B; dashed line), 323 kms" 1 (Model F15B; dash-dotted line), and 457kms _1 (Model H15B; dash-triple- 
dotted line). 
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Fig. 8. — Core hydrogen burning evolution in the HR diagram of 20 Mq models with ZAMS rotational 
velocities of Okms' 1 (Model D20; solid line), 103 kms" 1 (Model G20B; dotted line), 201 kms^ 1 (Model 
E20B; dashed line), SOTkms- 1 (Model F20B; dash-dotted line), and 425 kms" 1 (Model H20B; dash-triple- 
dotted line). The fastest rotating star becomes a Wolf-Rayet star at a central hydrogen abundance of 4%; 
the track of this star is discontinued at that point of evolution. 
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Fig. 9. — Evolution of the surface abundances during central hydrogen burning of 12 M Q (left column; 
Models D12, G12B, E12B, F12B, and H12B), 15 M (middle column; Models D15, G15B, E15B, F15B, and 
H15B), and 20 M (right column; Models D20, G20B, E20B, F20B, and H20B) stars for different ZAMS 
equatorial rotational velocities (see legend) as a function of time. Displayed are 12 C (Panels al-a3), 14 N 
(Panels bl-b3), and 16 (Panels cl—c3) relative to their initial abundance, and the 4 He mass fraction 
(Panels dl-d3). 
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Fig. 10. — Same models and abscissa as in Fig. 9, but the ratio of carbon to nitrogen (Panels al-a3), the 
ratio of carbon to oxygen (Panels bl-b3), and the ratio of nitrogen to oxygen (Panels cl-c3) are displayed, 
all relative to the abundance ratios in the initial models. In Panels dl-d3 the evolution of the equatorial 
rotational velocity is displayed and in Panels el-e3 that of the luminosity. 
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Fig. 11. — Same models and abscissa as in Fig. 9, but the evolution of the surface abundances of boron 
(Panels al-a3), 19 F (Panels bl-b3), and 23 Na (Panels cl-c3) are given relative to their initial values; for 
the radioactive isotope 26 Al (Panels dl-d3) the absolute mass fraction is displayed. 
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Fig. 12. — Evolutionary tracks in the HR diagram of rotating models in the mass range 8 Mq to 25 Mq (thick 
lines) and non-rotating 10 Mq to 25 M models (thin lines) from the ZAMS to the presupernova stage. The 
tracks are labeled with the initial masses (in M ). The tracks are for Models D10, D12, D15, D20, and 
D25 (no rotation) and for Models E08, E10, E12, E15, E20, E25 (ZAMS equatorial rotational velocities of 
~ 200 km s -1 ). The track of the rotating Model E25 is only shown until it evolves into a Wolf-Rayet star, 
during core helium burning. See also Heger et al. (1998). 
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Fig. 13. — Evolutionary tracks in the HR diagram of rotating models in the mass range 12 Mq to 20 Mq 
(thick lines) and of non-rotating models in the mass range 10 Mq to 25 M Q (thin lines) from core hydrogen 
ignition to the presupernova stage. The tracks are labeled with the initial masses (in M Q ). The tracks are 
for Models D10, D12, D15, D20, and D25 (non-rotating) and for Models E12B, E15B, and E20B (equatorial 
rotational velocities of ~ 200 kms -1 ). See also Fig. 12. 
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Fig. 14. — Evolutionary tracks in the HR diagram of stars with an initial mass of 12 Mq and ZAMS equatorial 
rotational velocities of Okms" 1 (Model D12, solid line), 99kms _1 (Model G12B, dotted line), 206 kms" 1 
(Model E12B, dashed line), 328 kms" 1 (Model F12B, dash-dotted line), and 474 kms^ 1 (Model H12B, dash- 
triple-dotted line). The evolution during central hydrogen burning is magnified in Fig. 66. The tracks start 
at the ZAMS and are followed until the phase listed in (Heger et al. 2000). 
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Fig. 15. — Evolutionary tracks in the HR diagram of stars with an initial mass of 15 M Q and ZAMS 
equatorial rotational velocities of Okms" 1 (Model D15, solid line), 102 kms -1 (Model G15B, dotted line), 
206 kms -1 (Model E15B, dashed line), and 323 kms -1 (Model F15B, dash-dotted line) up to core collapse. 
The evolution during central hydrogen burning is magnified in Fig. 7. 
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Fig. 16. — Evolutionary tracks in the HR diagram of stars with an initial mass of 15 M Q and ZAMS 
equatorial rotational velocities of Okms" 1 (Model D20, solid line), 103 kms -1 (Model G20B, dotted line), 
201 kms" 1 (Model E20B, dashed line), and SOTkms" 1 (Model F20B, dash-dotted line) up to core collapse. 
The evolution during central hydrogen burning is magnified in Fig. 8. 
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Fig. 17. — Surface mass fractions of 12 C (Panel a), 14 N (Panel b), and 16 (Panel c), relative to the 
initial values, as a function the initial equatorial rotation rate. Panel d gives the 4 He mass fraction. Solid, 
dotted and dashed lines give the values for the 12 M Q (Models D12, G12B, E12B, F12B, and H12B), 15 M Q 
(Models D15, G15B, E15B, F15B, and H12B), and 20 M (Models D20, G20B, E20B, F20B, and H20B) 
stars, respectively. The thin lines show the surface composition at the end of central hydrogen burning 
(Fig. 9), the thick lines those after the first dredge- up, i.e., when the star has first become a red supergiant 
after termination of central hydrogen burning. For the fastest rotating 15 M Q star we have no post-dredge-up 
data. 
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Fig. 18. — Equatorial surface rotation velocity (Panel a) and effective temperature (Panel b) for Model E12B 
as a function of time, displayed from the end of central hydrogen burning until core collapse. The dashed line 
indicates the effective temperature of the sun, and dotted lines separate "blue" from "red" phases. Between 
an age of 15.7 and 16.3 Myr the star undergoes a blue loop, while it spends the major part of central helium 
burning and the time from central helium depletion till core collapse as a RSG. Note that the maximum 
and minimum rotational velocity during the blue loop deviate by about one order of magnitude and that 
the rotational velocity after the blue loop is notably slower than before the loop. 
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